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ABSTRACT We introduce a Markov model for the gating of membrane channels. The 
model features a possible solution to the so-called gating current paradox, namely that 
the bell-shaped curve that describes the voltage dependence of the kinetics is broader 
than expected from, and shifted relative to, the sigmoidal curve that describes the voltage 
dependence of the activation. The model also predicts some temperature dependence of 
this shift, but presence of the latter has not been tested experimentally so far. 

INTRODUCTION 

The gating of membrane channels is of vital importance for the electrophysiological activity 
of nerve, heart and muscle. While some of these channels appear to have fractal-like gating 
(Liebovitch, 1995), most membrane channels do display activity that can be well approx- 
imated by a simple Markov process (Korn and Horn, 1988). However, Clay et al. (1995) 
revealed a gating current paradox that has been difficult to explain with a standard type 
(Hille, 1992) Markov model. The paradox is that the bell-shaped curve that describes the 
voltage dependence of the kinetics is shifted significantly relative to the sigmoidal curve that 
describes the voltage dependence of the activation. The standard type models (Hille, 1992) 
does not allow such a shift. Also the former curve is broader than the one predicted by the 
standard model. 

Here we introduce a new Markov model, that extends and generalizes the standard one. 
Our generalization consists of introducing an alternative route between the open and the 
closed positions of the gate. With two routes, or two membrane protein folding pathways, we 
are able to obtain results consistent with the observed ones. Thus such a model presents a 
possible resolution of the above paradox. A more complete resolution requires investigation 
of the detailed physical mechanism present in real membrane channels to see how they 
compare with the model. The idea with two routes, a rapid one and a slow one, is that the 
probability of choosing one or the other also depends upon the voltage through a Boltzmann 
factor. This will affect the kinetics, but not the equilibrium distribution (stationary state), 
and a relative shift of curves can take place. 

THE MODEL 

We imagine that a membrane channel has one open and one closed state as in the simplest 
standard (Hille, 1992) Markov model for this problem. However, between these states we 
now assume that there exist two routes (i = 1,2). This gives, 
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where the rate constants a±, a 2 and j3i, (3 2 , which are functions of voltage (but are constant 
at any given voltage), control the transitions between the closed (C) and the open (O) states 
of the gate. The ccj is the rate for a closed channel to open, and $ the rate for an open 
channel to close. We introduce effective rate constants a and (3, 

a = p l a 1 +p 2 a 2 (2) 
P = PiPi+P202 , (3) 

where the probabilities p 1 and p 2 are related in a standard way to the difference AGb in 
energy barriers that must be overcome for each of the two routes, 
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Let x denote the average fraction of gates that are open, or, equivalently, the probability 
that a given gate will be open, and let us imagine that a Markov (1906) model is suitable 
to describe the gating. One then has, as usual 

^ = a(l-x)-Px=— — -, (6) 
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Here Xoq denotes the steady stationary state fraction of open gates and r the relaxation 
time. At equilibrium, the probability for a channel to be in the open state is Xqo, and the 
probability to be in the closed state is (1 — x^). The ratio of these two probabilities is given 
by the Boltzmann distribution, 

ex P > ( 9 ) 
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where T is the absolute temperature, k is Boltzmann's constant, and AG X denote the energy 
difference between the open and the closed positions. Thus, 
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At equilibrium, each of the the forward reactions must occur just as frequently as each of 
the reverse reactions, giving, 
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This is the principle of detailed balance which is present in dynamical systems (reversible 
mechanics). As in the standard model the rates are then assumed to be, 

Q, ' = A - exp (w) (12) 

A = *«p(-1E5f), < 13 > 

where Aj is assumed to be independent of AG X . Thus the relaxation time (Eq. [8]) can then 
be written as, 
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a + P Picti + P2®2 + P1P1 + P2P2 
Using Eqs. g, §, [T2| and [T3|, we obtain 

2cosh(|§) 



(14) 



T = V 2kT I (Vo) 

(«! + A) exp(-|gt) + ( a2 + /3 2 ) exp(§§) 



cosh(f§) 



cosh(^) [A ieX p(-^) + A 2 exp(M|] 

cosh(^) 

A cosh(f§0 [exp(-fg| - 7) + exp(fg| + 7 
cosh(ff) 



2Acosh(f§0cosh(f§ + 7 ) 



(16) 
(17) 
(18) 



where, 
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To be more specific the voltage dependences of AG X and AG b are needed. For the energy 
difference between the open state and the closed state we assume as usual, 

AG X = Gciosod — G open = q x {v — v x ) — s x T , (21) 

where the term q x v x is due to the difference in mechanical conformation energy between 
the two states; q x v represents the electrical potential energy change associated with the 
redistribution of charge during the transition, and s x is due to the difference in entropy 
between the two states. A similar expression can be assumed for the energy difference 
between the two barriers in route 1 and 2, 

AG b = G 1 -G 2 = q b (v - v b ) - s b T . (22) 

Here v is voltage, while q x , v x , s x , q b , v b , and s b are constants. One notes that the curve for 
the relaxation time r has a shift in position due to the term 7. Inserted for the special case 
AG b = AG X the above yields, 



3 



./ x = I 1 + exp 
t = { 2X cosh 



(23) 
(24) 



-q x (v - v x ) + s x T 
kT 



-q x (v - Vx) + T(s x - 2-fk) 
2kT 



Here we find that the voltage dependence of the curve for the relaxation time (Eq. 
is shifted by an amount 2 r ykT/q x relative to the steady state activation curve (Eq. 
which means that the magnitude of the shift depends upon temperature. With AGb ^ AG X 
expression (Eq. |2"4] ) becomes more complex as follows from (Eq. |T8), and the shape of the 
former curve is modified. This however, is dealt with in the next section. 

RESULTS 

We now will compare the model with the experimental results of Clay et al. (1995) and show 
that it is consistent with the latter. Thus it presents a mechanism that represents a possible 
solution to the gating current paradox. The temperature-dependence of the currents were 
not considered in those experiments, so here s x and s& can be incorporated into v x and i>&. 
With use of Eqs. pll and |22|, Eqs. [TO and [18] become, 



1 



T 



1 + exp (=-=2 



cosh(^ 



2k b 



2Acosh(^)cosh(^ + 7 ) 



(25) 
(26) 



where k x = kT/q x and kj, = kT/q b . These expressions were evaluated numerically adjusting 
the parameters present to obtain a best possible fit to the experimental data. A least 
squares fit weighting various points in accordance with experimental uncertainty was used. 
The results of this evaluation is shown in the figure below where the data of Clay et al. 
(1995) is presented together with the curves given by Eqs. |25| and |26] using the parameters 
shown in the figure text. 

However, the curves are not very sensitive to the values of these parameters except 7, so they 
can be varied quite a bit and still give essentially the same curves. From these curves we 
find that the model is fully consistent with the experimental results within the uncertainties 
in the latter. 

DISCUSSION 

We have presented a Markov model that yields a possible solution to the gating current 
paradox announced by Clay et al. (1995). It gives a simple explanation of the voltage- 
shift of the bell-shaped curve for the relaxation time relative to the steady state activation 
curve. Also the width and shape of the relaxation time curve can be modified in a way 
consistent with experiments. A novel feature of the present model is that the voltage-shift 
is temperature dependent. It is not clear whether such a temperature dependence can be 
observed experimentally. 
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and 7 = 1.89. The error bars indicate the ± mean standard deviation from the six experiments of Clay et 
al. (1995). 
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